Absence of Two-Dimensional Bragg Glasses 
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The stability to dislocations of the elastic phase, or "Bragg glass" , of a randomly pinned elastic 
medium in two dimensions is studied using the minimum-cost-flow algorithm for a disordered fully- 
packed loop model. The elastic phase is found to be unstable to dislocations due to the quenched 
disorder. The energetics of dislocations are discussed within the framework of renormalization group 
predictions as well as in terms of a domain wall picture. 
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Randomly pinned elastic media are used to model var- 
ious condensed-matter systems with quenched disorder, 
including flux-line arrays in dirty type-II superconductors 
[1| and charge density waves [^J. Although it is known 
[3| that these systems cannot exhibit long-range transla- 
tional order in less than four dimensions, the intriguing 
possibility of a "topologically ordered" low-temperature 
phase remains an open question |^ ^| . It has been conjec- 
tured that such a phase, with all dislocation loops bound, 
exists in three dimensions. Such a phase would be elastic 
and have power-law Bragg-like singularities in its struc- 
ture factor; it is often referred to as a "Bragg" or "elastic" 
glass |§. 

Whether or not unbound topological defects exist 
at low temperatures involves a subtle balance between 
elastic-energy cost and disorder-energy gain ||. In this 
letter, we analyze this issue for two-dimensional ran- 
domly pinned elastic media at zero temperature by con- 
sidering a 2d-lattice model, viz., a fully-packed loop 
(FPL) model |jj with quenched disorder. Exact ground 
states of the FPL model are computed with and with- 
out topological defects, which we refer to as dislocations. 
The polynomial optimization algorithm jg] that we use, 
minimum-cost-flow, enables us to study large systems. 

We focus on the energetics of a single dislocation pair in 
systems of size LxL. Our main conclusion is that the dis- 
order energy gain of the optimally placed pair dominates 
over the elastic energy cost with the results being con- 
sistent with theoretical predictions of 0(— In 3 / 2 L) and 
0(+\tlL) respectively for the two quantities. Disloca- 
tions therefore become unbound and proliferate causing 
the destruction of the Bragg glass. 

Model and algorithm The FPL model is defined on 
a honeycomb lattice on which all configurations of oc- 
cupied bonds which form closed loops and cover every 
site exactly once are allowed, as shown in Fig. 1(a). This 
model can be mapped to a solid-on-solid surface model. 
Define integer heights at the centers of hexagons then ori- 
ent all bonds of the resulting triangular lattice such that 
elementary triangles pointing upward are circled clock- 
wise; assign 1 to the difference of neighboring heights 
along the oriented bonds if a loop is crossed and —2 oth- 



erwise. This yields single-valued heights up to an overall 
constant. It can be seen that the smallest "step" of the 
surface is three, so that the effective potential on the 
surface is periodic in heights modulo 3. 




FIG. 1. The fully packed loop model with periodic bound- 
ary conditions imposed. The ground states with and with- 
out a pair of dislocations for one realization of random bond 
weights are displayed in (b) and (a) respectively. The dislo- 
cations (solid dots) in (b) are connected by an open string 
(thick line) among the loops. The relevant physical object is, 
however, the domain wall which is induced by the dislocations 
as shown in (c) . This domain wall represents the line of bond 
differences between the ground states (a) and (b). 

Quenched disorder is introduced via random bond 
weights on the honeycomb lattice, chosen independently 
and uniformly from integers in the interval [-500, 500]. 
The energy is the sum of the bond weights along all loops 
and strings. The system can be viewed as a surface in a 
3d random medium that is periodic in the height direc- 
tion. Dislocations are added to the FPL model by "vi- 
olating" the constraint. One dislocation pair is an open 
string in an otherwise fully-packed system as shown in 
Fig. 1(b). The height change along any path encircling 
one end of the string is the Burgers charge ±3 of a dis- 
location so that the heights become multi-valued. Note 
that the configurations with and without a dislocation 
pair only differ on a domain "wall" , as shown in Fig. 1(c). 

Finding ground states is equivalent to an integer 
minimum-cost flow problem on a suitably designed 
graph; details can be found elsewhere (8||] . 

Numerical results For each disorder realization, the 
ground-state energies with and without a dislocation pair 
and hence the defect energy, Ed, and also the connecting 
wall (See Fig. 1(c)) were computed. 
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FIG. 2. Energetics of fixed defects: The probability distri- 
butions of the energy of a fixed dislocation pair with separa- 
tion L/2 for sample sizes from L — 12 to L — 384 are shown in 
(a). The solid lines are gaussian fits. The corresponding aver- 
age defect energy Ed (solid circle) and the root-mean-square 
(rms) width a(Ed) (solid square) are found to scale with sys- 
tem size as InL as shown in (b). The solid lines are linear 
fits. 

We first held the dislocations fixed at two specific 
sites separated by L/2 in L x L samples with L = 
12,24,48,96,192, and 384, with at least 10 4 disorder 
realizations for each size. The probability distributions 
P(Ed) of the defect energy are shown in Fig. 2(a); they 
are well fit by gaussians for all sizes. The mean pair- 
energy is plotted versus InL with a linear fit Ed = 
180(2) InL — 20(6). The root-mean-square (rms) width 
a (Ed) is also shown in Fig. 2(b), with a linear fit of 
a(Ed) = 250(3) + 133(l)lnL. This implies a tail in P(E d ) 
for negative Ed and suggests an energy gain from dislo- 
cations that can be optimized to take advantage of the 
negative part of the distribution. 

We have computed the optimized (lowest energy) dis- 
location pair energy E™™ for L up to 480 with 10 4 - 10 6 
samples for each size. The defect energy distribution 
P(E d uln ) is no longer Gaussian, indeed substantial asym- 
metry in P(E' d nin ) is seen in Fig. 3(a). Moreover, in con- 
trast to the case of fixed dislocations, E™ ln is negative 



and decreases more rapidly than InL with increasing 
system size while the rms width cr(L™ ln ) increases less 
rapidly than InL. The linear fits shown in Fig. 3(b) yield 
[\ E 7 [n \} 2/3 = (43-80 ± 0.31) ln(L) + (24.03 ± 0.14) and 
[cr(£™ in )] 2 = (21883 ± 180)ln(L) + (9013 ± 802). This 
implies that since almost all large systems have negative 
energy dislocation pairs, as is evident in Fig. 3(a) the FPL 
model is unstable against the spontaneous appearance of 
dislocations. 
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FIG. 3. Energetics of optimized defects: The defect en- 
ergy probability distributions for sample sizes from L = 12 to 
L — 480 are shown in (a). The solid lines are guides to the 

12/3 



eyes. Both the average defect energy plotted as 
vs. InL and the rms width plotted as [a(_E™ ln 
shown in (b). Solid lines are linear fits. 
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Continuum models To shed some light on the ob- 
served defect energetics, and to make contact with ana- 
lytic work, we consider coarse-grained continuum approx- 
imations to the random FPL model. In the absence of 
dislocations, the surface has a stiffness caused by the in- 
ability of a tilted and hence more highly constrained sur- 
face to take as much advantage of the low weight bonds 
as a flatter surface. In addition, the random bonds cou- 
ple to Vh as well as to h modulo 6 = 3. An appropriate 
effective Hamiltonian is thus 
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with f = f(r) locally random with variance A, and 
7 = 7(r) =0, ± 2 j- with equal probability. 

A related and well studied model (CDW) @ has, in- 
stead, the {7(1")} being uniformly distributed random 
variables in [0, 27r] . While it is not clear a priori that these 
models are in the same universality class, both the follow- 
ing, and an RG analysis, imply they are: shifting h(r) by 
^r7(r) changes the {/(r)} introducing short range cor- 
relations that are different in the two models, but these 
should, be irrelevant on large scales of h and L. The 
CDW model has an elastic glass phase for T <T g below 
which w is relevant and renormalizes to a T-dependent 
fixed point while A grows as InL yielding height varia- 
tions < [h{v) - h{0)} 2 > « ^Tln 2 r @. 

An approximate functional RG analysis yields a similar 
structure at all temperatures with a universal T = limit 
T of the coefficient T(T) [jjl|| . On large scales, the be- 
havior is dominated by the competition between the ran- 
dom stressing f(r), and the stiffness K, with the random 
force correlations effectively /j(q)/j(— q) = — Crying 
for small wave vectors q. We can thus work with the 
simple purely random- force limit of Eq. ([l]) with w = 0. 
We then immediately conclude that To = Ji b z ■ 

In the presence of dislocations h(r) becomes multi- 
valued. It can be decomposed into two parts h = h e + h d , 
with h e , the smooth elastic distortion while the singu- 
lar function h d , has a cut connecting the two disloca- 
tions at i"i and r 2 with Burgers charge b\ = b = 3 and 
b 2 = —b = -3 with V x Vh d — 2 W 2 (r — r.;) and 

V 2 h d = 0. The energy of a dislocation pair is 



E d = ^- I d 2 r\Vh, 12 



2 
Kb 2 
~2^~ 



J d 2 r(f T ■ Vh d ) (2) 
Inlrt-ral-ftfertrO-stra)] . (3) 



where the static random force field has been decom- 
posed into longitudinal i L = Vu(r) and transverse f T = 
V x g(r) components. Since u(r), g(r), and Vh d are con- 
tinuous across the cut while h d jumps by b, f L makes no 
contribution to E d and Eq.(^) follows by integration by 
parts. The first term is the standard elastic cost and the 
second the disorder gain so that g(r) is the potential felt 
by the dislocations. Its variance is S g = p(q)<?(— q) = 
— Cq~ 2 \nq 2 so that the elastic surface without disloca- 
tions and the dislocation potential have the same statis- 
tics after a rescaling: Sh = h(q)h(—q) ~ S g /K 2 . 

In terms of the statistical properties of the dislocation 
potential our numerical results can be understood. We 
first discuss fixed dislocation pairs. It has been argued 
jnj from a statistical symmetry (of the CDW model) that 
on large scales /(r) and hence the potential g(r) and E d 
are gaussian. The shape of our computed defect energy 



distribution, its mean and its variance all agree with pre- 
dictions from Eq. (||), E d w (K b 2 /2ir) InL — an exact 
result for the CDW model— and a{E d ) « y/WU/nlaL 
with K = 126(2) and C = 6174(93) yielding C/K 2 = 
0.389(11). We also measured the variance of the height 
of the surface without dislocations finding a good fit 
to a 2 (h) = 0.061(2) In 2 L + 0.477(7) InL + 0.765(13). 
The predicted coefficient of the In 2 L is C/2irK 2 yield- 
ing C/K 2 — 0.387(12). The agreement between the two 
estimates of C/K 2 further supports the validity of the 
random force model, including, in particular, the equal- 
ity between the longitudinal and transverse parts of f (r) 

We now turn to optimized dislocations. If the first 
term in Eq. @ can be ignored, the energy of a dislocation 
pair will be lowest if the dislocations are at the minimum 
(g m i n ) and maximum (g max ) of the random potential, 
g{r). The distribution of the extrema of potentials like 
g(r) whose variance grows as a power of InL, can be 
semi-quantitatively understood by thinking of iteratively 
optimizing over each factor of two in length scale, with 
the component of the random potential on scale I be- 
ing essentially the contribution from Fourier components 



with 



^ < \Qx\,\q y 



< 2v ? 7r . If scale I gives rise to a 



contribution to the variance of g of order (In l) 2a — with 
a = i in our case — , then a typical g(r) is the incoherent 
logarithmic sum over all scales, i.e. of order (In l) a+1/>2 . 
The maximum of g can be found, heuristically, by maxi- 
mizing over the four points at scale I = 1 in each square 
of scale 2, then maximizing over the four scale 2 maxima 
in each scale 4 square, etc. Since the scale I structure 
of g is weakly correlated over scales much longer than 
I, a crude approximation is to ignore these correlations 
whereby one stage of optimization adds of order (In l) a 
to the local scale I maximum [^2|. Thus scales should 
be summed over coherently yielding g max ~ (lnL) a+1 . 
From this hierarchical construction, it can be seen that 
the variance of <7 max is dominated by the largest scales, 
so that cr(g max ) ~ (lnL) Q . In our case, we thus expect 
ffmax — ^min ~ (InL) 3 / 2 which indeed dominates over the 
InL elastic energy term in Eq. (|J). Hence we expect 
Ef n w -AbVC ; (ln L) 3 / 2 and cr(L^ nin ) w BbVCilnL) 1 / 2 
with some coefficients A and B. 



The linear fits in Fig.3(b) give A » 1.23(1) and B w 
0.56(1) using the C computed earlier. If the elastic part 
In L is subtracted from L™ ln by fitting the difference 
L™ in (L) - L™ in (L/2) to - (31n2/2)Afe v / C(lnL) 1 / 2 + 
const., this yields a very comparable A w 1.17(5). Ex- 
tremal heights hf in = h min — h max of the elastic sur- 
face without dislocations and optimal dislocation ener- 
gies can also be used to extract K via L™ ln — bKh™ ln w 
Tp- InL + const.. This yields K ps 114(1) in not unrea- 
sonable agreement with the K s» 126(2) from the fixed 
dislocation pairs. 

An upper bound on gmax — 3min can be simply ob- 
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tained by noting that Prob[g max > M] < Prob[g(r) > 
M] so that, with L 2 points, the median <? max is less than 
the M for which the right hand side, L 2 Prob[g(r) > M] 
for fixed r, equals \. If g(r) is gaussian sufficiently far 
into the tail of its distribution, this yields g max — 3mm < 

VC x f^(lnL)i so that A < \f^- = 1.5957.... The hi- 
crarchical optimization described above suggests that A 
might saturate this bound. To test the universality of 
A and this issue, we have measured the distributions of 
the extrema of the heights of the FPL surface without 
dislocations and several simulated exactly gaussian ran- 
dom surfaces with different a. Good fits are found to 
(In L/a a ) 1+a with a a a cutoff, yielding Afpl heights ~ 
1.45(3), and A gaussian « 1.517(3), 1.307(6), 1.168(6), and 
1.064(5) for a — 0, 1/2, 1, and 3/2 respectively. Similarly 
from the variances we obtain -Bfpl heights ~ 0.67(2), and 
B gaU ssia„ Pa 0.475(4), 0.637(7), 0.730(9), and 0.814(4) for 
a = 0, 1/2, 1, and 3/2 respectively. 

We first note that all the extracted A's satisfy A < 
and there appears to be a systematic trend for 
A(a) to decrease as a increases; if this is true then it is 
most likely that A is strictly less than the bound for all 
a > — i (for a < — i gaussian surfaces the variance sat- 
urates for large L and the extrema grow as VhiL). The 
values of A for the nominally a = ^ cases, 1.17-1.23, 1.45 
and 1.307 for the optimal dislocations, extrema heights, 
and the gaussian surface respectively, differ by substan- 
tially more than the apparent statistical errors as do the 
B's, 0.56, 0.67, 0.637 respectively. But given the narrow 
range available of In L in spite of a large range of sizes 
and, as importantly, the lack of understanding of cor- 
rections to scaling, these results are certainly consistent 
with universal values of A and B for a = \. At this 
point, however, understanding whether this is in fact the 
case, and also whether A for gaussian surfaces depends 
on a, must wait for better theoretical understanding. 

Overall, we have found rather good agreement for a 
variety of large scale quantities with the RG prediction 
of equivalence at long scales of the FPL model and a 
random force model. Although extracting reliable expo- 
nents of In L is not possible (especially with logarithmic 
corrections to scaling) the fact that the coefficients and 
ratios between these — extracted several ways — are in rea- 
sonable agreement is a more stringent test. But even if 
the random force equivalence is not valid, the data of 
Fig. 3 clearly indicate the instability of large systems to 
dislocation pairs. With no restrictions on their number, 
dislocations will proliferate thereby driving the elastic 
constant K to zero. 

We conclude with an alternate way to understand the 
structure of excitations in the elastic glass, via a pic- 
ture developed for the three-dimensional case Q. The 
basic excitations from a ground state are fractal domain 
walls surrounding regions in which h changes by b. Their 
fractal dimension, d w , will be the same as that for the 



forced open wall that connects a pair of dislocations 
(Fig. 1(c)) for which we find d w — 1.28(3) for fixed pairs 
and 1.30(3) for optimized pairs. [These contrast strongly 
with the connecting strings in the loop model which have 
d s = 1.75(3) and 1.74(3) respectively, very close to the 
value in the non-random loop model 0]. The energy of 
a scale L wall constrained only on scale L is predicted 
to vary by of order Vhi L but have mean independent 
of L. The incoherent logarithmic addition over all scales 
then yields variations of the fixed-end open domain wall 
energy, of order In L and a mean of the same order — as 
found. But if the end positions can adjust to lower the 
wall energy near the dislocation at each scale, the en- 
ergies add up coherently resulting in the — ln 3//2 L mean 
optimal dislocation pair energy with order y/ln L around 
the mean variations being dominated by the largest scale, 
in an analogous way to the extrema of the random po- 
tential g{r) of the random force picture. Since the defect 
energy in the domain wall picture is concentrated on the 
wall, while it is spread out over a region of area 0(L 2 ) in 
the random force model, it is surprising that these yield 
the same predictions! But the fact that our results agree 
well with the domain wall picture in 2D lends strong sup- 
port to the validity of the analogous picture in the 3D 
case for which it has been used to conclude that the 3D 
elastic glass phase is stable to dislocation loops ||. 
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